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FOREWORD 


The present tables were prepared by the Applied Mathematics 
Research Laboratory of Aerospace Research Laboratories, under 
Project No.7071. ‘The computations were initiated by Mr. Henry E. 
Fettis, and carried out by Mr. James C. Caslin on the IBM 1620 


and 7094 electronic digital computers. 


ABSTRACT 
This report contains 15 place tables of the modified 


Bessel Functions I (x), I, (x), eT (x), e*I, (x), for 


x = 0(.001)10. 
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The function T(z) is 


equation 


dy 


2 
dy , L 
z2 az 


a 


which is finite at z = 


It has the following power 


form 


INTRODUCTION 


that solution of the modified Bessei 


ye 
yt oy 0 
z 


(1) 


O and which behaves in this neighborhood like 


(2) 


series development 


z 2 z,4 
5)” ) ) 
1 ee l+ ———— + +...4 (3) 
P(v + 1) (v + 1)-1! (vi + 1l)-(v + 2)-2! 

and for v =n, an integer, the integral representation 

I (2) = = fre © 808 © cog: nb: a0 (4) 
When lz| is large, the function T(z) behaves like 

ry Zz 
lier e (5) 


and for this reason it is sometimes preferarle to tabulate it in the 


I (2) 


METHOD GF COMPUTATION 


The tables were computed from the integral representation 
{equation 4), making use of the fact that if the integral is 
approximated by the trapezoidal rule using (m+ 1) ordinates, 


the error is of the order of magnitude of 


Table (2) lists values of the quality I (2) for selected values 

of z and v. Inspection of this table enables one to predict the 
accuracy attainable for a given number of ordinates. In the 

present compi?ation a maximum of 32 ordinates was required to give 
better than 16 places of accuracy in e* I (z) and e~ I, (z) for 

z< 10. As z became smaller, a lesser number of ordinates was 
needed. The nuinber 32 was chosen since it enabled the same ordinates 
to be used throughout by successively having the interval of integra- 
tion with increasing z, starting with only two subdivisions for the 


smaller values. 
COMPARISON WITH OTHER TABULATIONS 


The most accurate computations of the modified functions were 
done sometime ago by Aldis [1]. dis results have proved surprisingly 
accurate, particularly where one considers the limited tools available 


at that time. Comparison of the present values with Aldis' indicates 


that the last digit is inaccurate in some cases. This is due to 

the fact that a maximum of 16 digits were available on the computer 

with the unavoidable introduction of roundoff errors (see [2]). 
Comparison of the present 16 place values with the 10 place 

ones in [3] have disclosed numerous inaccuracies in the final digits 

of the latter tabulation. These have been listed in a separate note 


[4] and are repeated here in Table 3 for the sake of completeness. 


(1) 


(2) 


(3) 


(4) 
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